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Abstract. The electronic structures of the metallic and insulating phases of the alloy series Ca2- a; Sr 2 ,Ru04 
(0 < x < 2) are calculated using LDA, LDA+U and Dynamical Mean-Field Approximation methods. In 
the end members the groundstate respectively is an orbitally non-degenerate antiferromagnetic insulator 
(x = 0) and a good metal (x — 2). For x > 0.5 the observed Curie- Weiss paramagnetic metallic state which 
possesses a local moment with the unexpected spin S — 1/2, is explained by the coexistence of localized 
and itinerant Ru-4d-orbitals. For 0.2 < x < 0.5 we propose a state with partial orbital and spin ordering. 
An effective model for the localized orbital and spin degrees of freedom is discussed. The metal-insulator 
transition at x — 0.2 is attributed to a switch in the orbital occupation associated with a structural change 
of the crystal. 

PACS. 7 5.30.-m, 75.50.-y, 71.25.Tn 



1 Introduction 

The discovery of unconventional superconductivity in 
Sr2RuC>4 |y|2J has evoked considerable interest in the elec- 
tronic properties of ruthenates. Curiously the substitution 
of the smaller Ca 2+ -ions for Sr 2+ -ions does not lead to 
a more metallic state but to an antiferromagnetic (AF) 
Mott insulator with a staggered moment of S = 1 as ex- 
pected for a localized Ru 4+ -ion which has 4 electrons in 
the £2g-subshell. As will be discussed below, this insulat- 
ing behavior is driven by a crystallographic distortion and 
a subsequent narrowing of the Ru-4d bands. The com- 
plete series of isoelectronic alloys for intermediate con- 
centrations has recently been synthesized and studied by 
Nakatsuji and Maeno [HQ]. This gives a rare opportunity 
to examine the evolution of the electronic structure from a 
multi-band metal to a Mott-insulator transition in an iso- 
electronic system. The evolution is not at all monotonic 
and the metal-insulator transition does not take place as 
a simple Mott transition but proceeds through a series of 
intermediate regions with unexpected behavior. 

The most dramatic example is the system at a con- 
centration x = x c = 0.5. At this critical concentration 
the susceptibility shows a free Curie form with aS = 1/2 
moment (not S = 1) per Ru-ion coexisting with metal- 
lic transport properties. This critical concentration repre- 
sents the boundary of the paramagnetic metallic region 
which evolves as Ca is substituted in the good metal, 
Sr 2 RuC>4. At higher Ca concentrations (x c > x > 0.2), the 
alloys enter a region with AF correlations at low tempera- 



ture but still with metallic properties. Insulating behavior 
appears only at smaller values of x < 0.2. The challenge 
that we address in this paper, is to understand this unex- 
pected and nonmonotonic evolution and in particular the 
exotic behavior in the vicinity of the critical concentration 

We begin by discussing the end members. Sr2Ru04 is 
a good metal. Its Fermi surface has been determined by de 
Haas- van Alphen experiments Jjj and agrees very well with 
the predictions of the local density approximation (LDA) 
in the density functional theory M . Ca2RuC>4 is a AF Mott 
insulator and it can be well described by augmenting the 
LDA by a mean field to include the onsite correlation — 
the so-called LDA+U method M . The third section of this 
paper is devoted to the examination of the intermediate 
concentrations by performing calculations for a series of 
characteristic x-values. The most challenging is the region 
x w x c = 0.5 where there are strong correlations but no 
symmetry breaking so that both LDA and LDA+U are in- 
applicable. At this concentration we employ the recently 
developed the ab initio computational scheme combining 
local density approximation and d yna mical mean field the- 
ory @,0Jl| (LDA+DMFT) ||||. We use LDA calcu- 
lations to determine the input parameters in the effective 
Anderson impurity model which in turn is treated using 
a non-crossing approximation (NC A) Ea] . These calcula- 
tional schemes give us reliable information on the evolu- 
tion of the electronic distribution among the 3 orbitals in 
the i2g-subshell, which we shall show is the key to under- 
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1. Basic crystall structure of isoelectronic alloy series 
E Sr 2 ,Ru04. 
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Fig. 2. Local electronic structure of isoelectronic alloy series 
Ca2-iSra;Ru04. In Ca2Ru04 spin-down electron occupies xy- 
orbital (left panel); In S^RuCU spin-down electron occupies 
a;,z/yz-orbitals (right panel); 



standing the electronic properties. The paper concludes 
with a discussion and summary of our results. A brief ac- 
count of this work has appeared elsewhere |i~6| . 



2 End Members: Sr 2 Ru0 4 and Ca 2 Ru0 4 

We start with Sr2Ru04 (or x = 2). This is a good metal, 
forming a 3-dimensional but anisotropic Landau- Fermi liq- 
uid at low temperatures OJl8[. Sr2Ru04 crystallizes in 
the undistorted single-layered K^NiF^structure U^Of (see 
Fig. [j]) with lattice parameters quoted in Table yTThe 
Ru06-octahedra are slightly elongated along the c-axis. 
The Ru-ions have a formal valence Ru 4+ and have a tetrag- 
onal local symmetry. The 2p-0 levels are completely filled, 
leaving 4 electrons in ^g-subshell of the 4d-Ru levels. The 
crystal field level scheme that would apply for an isolated 
Ru 4+ -ion is shown in Fig. g. The upper e s -shell (not in- 
cluded in this figure) is empty. The splitting between the 
xy-orbitals and the degenerate {xz, j/z}-orbitals is small. 
But the xy-orbitals 7r-hybridize with 2p-orbitals of all 4 in- 
plane O-neighbors while the xz(y2:)-orbitals 7r-hybridize 
only with the 2 O-neighbors along the x(y)-axis. As a re- 
sult the xy-bandwidth is approximately twice the {xz, yz} 
bandwidth (see Fig. ||). The LDA calculations || give 3 
Fermi surface sheets, one with essentially xy and two with 
mixed {xz,yz} character. Their shape and volume agree 
with the de Haas- van Alphen results J7] . 

The volumes contained by the Fermi surface sheets 
give an almost equal occupancy of each of the 3 £ 2g -orbitals. 
If we denote the occupancy of the {xz,yz} and (xy)- 
orbitals by (nr a m,n-y), then Si'2Ru04 has the fractional 
occupancy (8/3, 4/3). Altho' there are clear signs of strong- 
correlations in the enhanced effective mass (enhancements 
~ 3 — 4 [[7 pl| , |22| ] ) and low effective Fermi temperature, the 
low-temperature behavior is clearly that of a well-defined 
Landau-Fermi liquid. 

Turning to the other end member, Ca2RuC>4 or x = 0, 
the substitution of the smaller Ca 2+ -ion for Sr 2+ does not 




(8/3, 4/3) 



Energy (eV) 

Fig. 3. Density of tig states for Sr2Ru04 obtained from LDA 
calculation. The solid line is the DOS for the (xz, j/z)-orbitals 
and the dashed line for the xy-oibital. (n( yZiZX -\,n X y) indicates 
the electron occupation of the orbitals. 



lead to a uniform shrinking of the lattice parameter. In- 
stead the RuOg-octahedra undergo a combined rotation 
and tilt (P&ca-structure) so that the Ru-0 bond length is 
preserved but the Ru-Ru separation contracts. In Fig. [|we 
illustrate the relevant distortion of the crystal structure. 
This distortion bends the Ru-O-Ru bond angle away from 
180°, thereby reducing the bandwidth of the t2g-orbitals. 
Also the smaller size of the Ca 2+ -ion decreases the in- 
terlayer distance (i.e. the c-axis lattice constant) which 
results in a change from elongation to a compression of 
the RuOg-octahedra. This in turn changes the sign of the 
energy splitting between the (xy)- and (xz, yz)-orbitals, so 
that now the xy-orbital lies lower in energy (see Fig. |2|). 
The crystal structure is orthorhombic (see Table [l]) . All 
RuOg-octahedra are equivalent with a rotation around 
their long axis (0 1) and a tilt around the diagonal in- 
plane axis (110) (Fig. ||d). Note all inplane O-ions are 
equivalent in this structure. 
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a) Sr 2 Ru0 4 



i)Ca ]5 Sr O gRii0 4 




c) Ca 18 Sr 02 RuO 4 



d) Ca 2 Ru0 4 



Fig. 4. Scheme of crystall distortion of Ca2- a: Sr a; Ru04 Con- 
secutive structural change of the O-octohedra in the alloy series 
Ca2- a; Sr I Ru04: a) Ideal structure K2NiF4-type (space group 
I4/mmm); b) Space group I4i /acd derives from 14/ mmm by 
rotation around [001]-axis; c) Space group P2i/c described by 
the additional rotation around a free axis in the octahedron ba- 
sis plane, d) Space group P&ca derived from the ideal structure 
by rotation around the [001]- and [110]-axes. 



Ca 2 Ru0 4 is an AF insulator. The LDA+U method @ 
which is based upon spin-orbital unrestricted Hartree-Fock 
equation (i.e. a static mean field treatment), generally 
works well for magnetic long range ordered insulators |23f . 
We applied this method to Ca2RuC>4 choosing parame- 
ter values U — 1.5 eV and J = 0.7 cV for the onsite 
Coulomb repulsion and the intra-atonric Hund's Rule cou- 
pling. The method converged to an AF insulating ground- 
state in which the lower icy-orbital is fully occupied and 
the majority spin (xz, yz)-orbitals are also occupied. The 
sublattice magnetization is reduced by hybridization with 
the O-orbitals from the full value of 2 hb expected for 
S — 1, to 1.35 /is. The energy gap is unusually small 
(0.17 eV). Both these values agree well with experiment, 
1.3 tis [[19) and 0.2 eV |gj| respectively. Further treating 
spin-orbit coupling in second order perturbation theory 
gives preferred orientation of the magnetic moment along 
the orthorhombic 6-axis (or [1 1 0] in tetragonal notation) 
also in agreement with experiment [£5| (see Appendix). 

The LDA+U method gives a satisfactory description 
of the electronic structure of Ca2RuC>4. In terms of our 
previous notation, Ca2RuC>4 has an integer orbital occu- 
pancy of (2,2). The key issue that we address below is the 
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Fig. 5. Phase diagram of Ca2_^Sr a; Ru04. The phases are de- 
noted by the abbreviations P-I for paramagnetic insulating, P- 
M for paramagnetic metallic, AF-I for antiferromagnetic insu- 
lating, AF-M for antiferromagnetic metallic, M-M for magnetic 
metallic and SC for superconducting (Figure is taken from M) 



evolution of the electronic structure between these very 
different end members. 



3 Evolution of the electronic structure at 
intermediate concentrations 

The results of our calculations presented above confirm 
the contrasting behavior of the end members, Sr2Ru04 
and Ca2RuC>4. The experimental investigations of Nakat- 
suji and Maeno show a complex evolution of the elec- 
tronic structure which they break down into 3 concen- 
tration regions, labelled Region I (0 < x < 0.2), Region II 
(0.2 < x < 0.5) and Region III (0.5 < x < 2), each with 
its own characteristic behavior. Their results are summa- 
rized in the phase diagram (Fig. ||). Below we discuss each 
of these regions, starting with the Sr-rich region. 



3.1 Region III (2 > x > 0.5) 

The experiments of Nakatsuji and Maeno show that super- 
conductivity is rapidly suppressed when Ca is substituted 
for Sr in Sr'2Ru04. This suppression is a natural conse- 
quence of disorder for an unconventional superconductor. 
More interesting is the evolution of the spin susceptibil- 
ity which progressively increases with Ca substitution and 
evolves from an enhanced Pauli susceptibility into a Curie- 
Weiss form. The characteristic (Curie- Weiss) temperature, 
0cw{< 0), approaches zero for x — > 0.5. At the same time 
the linear term in the low temperature specific heat is 
also enhanced but the Wilson ratio defined at low tem- 
peratures is strongly enhanced and approaches a value of 
40 as x — > x c . 

This value of x c = 0.5 is a critical value, separat- 
ing a metallic and orbitally ordered phase with antifer- 
romagnetic spin correlation (x < x c ) from the paramag- 
netic metal for x > x c . As mentioned above, one observes 
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9 CW rs for x s=s x c . Even more remarkable is the evolution 
of the magnetic moment which takes a value of S = 1/2 
as x — > x c (5|(| ■ This value is quite distinct from the value 
of S = 1 in Ca2RuC>4 — the expected value for a localized 
Ru 4+ -ion with 2 holes in the i2g-subshell. Moreover, alloys 
with x ~ x c are metallic, not insulating. The explanation 
of this unexpected behavior presents a clear challenge. 

As discussed above Sr 2 Ru04 is a good metal with an 
electronic structure that divides into two distinct and frac- 
tionally filled bands, a d^-band with f» 4/3 electrons per 
Ru and the hybridized d XZt y Z -bajid pair with 8/3 electrons 
per Ru. Now the evolution with Ca substitution towards 
a Mott insulating state occurs because the bands are nar- 
rowed by rotation of the RuC>6 octahedra. In general, the 
evolution from a metal to a Mott-insulator is driven by 
growing Umklapp scattering. This is very clear in one di- 
mension where studies of chains and ladders show insu- 
lating behavior at half-filling already for arbitrarily small 
Coulomb repulsion driven by elastic Umklapp scattering 
processes across the Fermi surface. Here we are dealing 
with an approximately two-dimensional multi-band situa- 
tion. In the single band case the Umklapp surface, across 
which elastic Umklapp processes are allowed, will in gen- 
eral be different from the Fermi surface, but at half-filling 
each encloses the same volume. In Sr2Ru04 the fractional 
occupancy of each subband clearly forbids elastic Umk- 
lapp scattering across the Fermi surface in low orders. In- 
deed this may well be the reason that Sr2Ru04 forms a 
good Landau-Fermi liquid, even though it is clearly close 
to a Mott insulator. One way to enhance Umklapp scat- 
tering is to transfer electrons between the subbands. This 
will allow low-order Umklapp scattering, when integer oc- 
cupancy of the subbands is reached. 

The theoretical investigation at these concentrations 
is inhibited by the lack of symmetry breaking at x ~ x c . 
Therefore we cannot use the LDA+U method to build in 
the onsite correlations that are totally neglected in stan- 
dard LDA calculations. Recently considerable progress has 
been achieved on the theory of the Mott transition in the 
Hubbard model by the use of the dynamic mean field the- 
ory (DMFT) method |lJ|L2|. This is essentially an ex- 
pansion around an infinite coordination number and for- 
mulates the problem in terms of an effective Anderson 
impurity model which is to be solved self-consistently. In 
this way the growth of onsite correlations can be treated 
as the Mott transition is approached in a paramagnetic 
metal. Recent advances use LDA calculations to determine 
the input parameters and a non-crossing approximation 
(NCA) to solve the effective Anderson model. 

We performed a series of calculations usingthis LDA 
+ DMFT (NCA) approximation scheme p6p7| for the 
Sr 2 RuC>4 structure. We increased the value of Hubbard- 
U to examine how the onsite correlations grow. In Fig. g 
is a series of results for the density of states (DOS) in 
the xy- and (xz,y z)-subbands. Since these subbands have 
quite different widths, the onset of Mott localization oc- 
curs at different critical values of U. Thus we see that as 
U is increased through a value of U « 1.5 eV there is 
a transfer of electrons between the subbands so that the 
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Fig. 6. Results of LDA+DMFT(NCA) calculations obtained 
within LDA DOS for Sr 2 Ru0 4 . The solid line is the DOS for 
xz, yz-orbitals and the dashed line for (xy)-OTbital. At U = 
1.5 eV the xz, yz-orbitals become localized. At U — 2.5 eV 
additionally the localization of sy-orbital occurs. The Fermi 
energy is defined to be zero and was adjusted to conserve the 
number of particles (4 electrons per site) . 



integer occupancy of 3 electrons and Mott localization ap- 
pears in (xz, yz)-subbands while the broader half-filled xy- 
band remains itinerant. This unusual behavior is driven by 
the combination of the crystal field splitting, as shown in 
Fig. @ [(xz, yz) lower] and the narrower bandwidth of the 
(xz, j/z)-orbitals. A further increase in the value of U to 
U rs 2.5 eV is required to obtain Mott localization also in 
the sy-subband. 

These results lead us naturally to the following pro- 
posal to explain the anomalous properties in the critical 
concentrations x = x c . The electronic configuration is now 
(3,1). The 3 electrons in the {xz,yz}-subbands are Mott 
localized and have a local moment of S = 1/2. The re- 
maining valence electrons are in the itinerant xy-ha,nd and 
is responsible for the metallic character. Thus at this con- 
centration we have the unusual situation of localization 
in only part of the 4d-orbitals and coexisting localized 
and itinerant 4d-orbitals. Note that in the orthorhombic 
crystal structure at x — x c the 2 subbands have differ- 
ent parity under reflection around a Ru02-plane, similar 
to tetragonal Sr2Ru04, which forbids direct hybridization 
between the subbands. This proposal explains in a natural 
way the unexpected moment of S = 1/2 of the Ru-ions and 
the coexistence of metallic behavior and local moments. 

Note that the calculations are carried out more con- 
veniently by increasing the value of the onsite repulsion, 
U which however should not change appreciably with the 
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Table 1. Crystallographic data, which is used for LDA calculations: symmetry group, parameters of lattice, atomic positions 
and distance between nearest atoms. Symbol "/" denotes that for this structure corresponding parameter do not exist. 



Compound 


Sr2RuOi 


Cai.sSro.5Ru04 


Cai.$SrQ.2RuC>4 


Co^RuOa 


Symmetry group 


14/mmm 


IAi/acd 


P2i/c 


Pbca 


o[A] 


3.8603 


5.3195(1) 


5.3338(4) 


5.6323(3) 


MA] 


3.8603 


5.3195(1) 


5.3162(4) 


11.7463(5) 


c[A] 


12.729 


25.1734(5) 


12.4143(8) 


5.3877(2) 


Voi.lA 3 } 


189.69 


712.33(2) 


352.01(4) 


356.45 


/?["] 


/ 


/ 


90.06(1) 


/ 


Ca(Sr) x 


0.0 


0.0 


0.0141 (21)/0.4903(24) 


0.0593(4) 


Ca(Sr) y 


0.0 


0.25 


0.0137(23)/0.5273(23) 


0.3525(2) 


Ca(Sr) z 


0.14684 


0.5492(1) 


0.3483(2) 


0.0021(5) 


Ol x 


0.0 


0.1933(2) 


0.1939(6) 


0.3005(4) 


Ol y 


0.0 


0.4433(2) 


0.3079(6) 


0.0272(2) 


Ol Z 


0.3381 


0.125 


0.0/0.0196(5) 


0.1952(4) 


02 X 


0.5 





-0.0344(5) 


-0.0212(4) 


02 y 


0.0 


0.25 


-0.0064(7) 


0.1645(2) 


02 z 


0.0 


0.4568(1) 


0.1649(2) 


-0.0692(3) 


Ru-01[A] 


1.930 


1.929(1) 


1.936(3)/1.926(3) 
1.941(3)/1.952(3) 


2.015(2) 
2.018(2) 


Ru - 02 [A] 


2.061 


2.059(3) 


2.056(3)/2.056(3) 


1.972(2) 


Co -01 [A] 


2.692 


2.399(2) 


2.316(7)/2.286(10) 


2.292(3) 






2.994(2) 


2.445(8)/2.502(9) 
2.838(11)/2.934(10) 
3.141(10)/3.037(10) 


2.433(3) 
2.565(3) 
3.313(3) 


Ca - 02 [A] 


2.439 


2.326(4) 


2.294(4)/2.296(4) 


2.287(3) 




2.737 


2.664(1) 


2.416(12)/2.488(13) 
2.559(13)/2.444(13) 
2.772(13)/2.845(13) 
2. 932(12)/2. 912(13) 


2.362(3) 
2.399(3) 
3.118(3) 
3.296(3) 



concentration, x. In reality it is the bandwidth which is 
changing with the decreasing x as the RuOg-octahedra 
progressively rotate when Ca is substituted for Sr. The 
key result however is the existence of a parameter range 
where this partial localization is stable. The fact that we 
calculated only for the highly symmetric Sr2Ru04 struc- 
ture, rather than the distorted structure is, we believe, 
unimportant in establishing this (3,1) configuration as a 
stable electronic configuration. 



3.2 Region II (0.5 > x > 0.2) 

At lower values of x we enter Region II (0.5 > x > 0.2) 
characterized by a tilting plus rotation of RuCvoctahedra. 
Ca1.8Sro.2R.uO4 has a low-symmetry crystal structure with 
the space group P2i/c [E8j, which can be obtained from 
the tetragonal IA/mmm structure by rotating and tilt- 
ing of the RuC>6-octahedra similar to pure Ca2Ru04 but 
with a smaller tilting angle [E8| (Fig. He). There are now 
two types of in-plane oxygen ions and two types inequiv- 
alent of RuCvoctahedra. The RuOg-octahedra continue 
to be elongated in this region so that the xy-orbital con- 
tinues to lie higher in energy. The metallic character of 
the alloys in this region shows that the itinerant character 
of the sy-subband is preserved, although the bandwidth 
will be narrowed by the additional tilting distortion of the 
RuOg-octahedra. Our conclusion is that the (3,1) orbital 



occupation continues to hold also in Region II with local- 
ization of the electrons only in the {xz, yz}-subband. 



3.3 Region I (0.2 > x > 0) Ca-rich 

The Ca-rich region is characterized by a transition to an 
insulating groundstate and simultaneously a change in the 
crystal structure. The S-Pbca structure of the groundstate 
in this region evolves continuously out the groundstate of 
pure Ca2Ru04. As discussed above, the change from elon- 
gated to compressed RuOe-octahedra causes a switch in 
LDA+U in the orbital occupation numbers to (2,2) i.e. 
to a filled xy-subband and a half-filled {xz, yz}-complcx. 
Therefore we assign the insulating groundstate in all of 
Region I, to an orbital occupation (2,2). The first order 
transition between Regions II and I as a: is decreased be- 
low 0.2, is to be associated with a switch in the orbital 
occupation from (3,1) to (2,2). 



4 Magnetic Properties of the different regions 

4.1 Region III 

We now turn to the magnetic properties at x > x c — 0.5. 
These are dominated by the local S — 1/2 moments. A 
Kondo-type of interaction between the two bands [xz, yz) 
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and xy can be excluded due to the absence of hybridiza- 
tion. The xy-orbital cannot hybridize with the other two, 
since they have opposite parity with respect to reflection 
on the basal plane. Via Hund's rule coupling, however an 
RKKY interaction between the localized spins is induced, 



Wrkky 



q 



j£ X (q)S q ■ S_ c 



(1) 



where S q = S yZi q + S zx ^ and x(q) is the static spin sus- 
ceptibility in the d^-band and Jh is the onsite Hund's 
rule coupling, and promotes antiferromagnetic correlations. 
At the same time, however, Hund's rule coupling may 
cause the ferromagnetic correlations through the double 
exchange mechanism. Thus, the two types of spin interac- 
tions mediated by the itinerant electrons of the a;y-band 
tend to compensate each other, such that the net ex- 
change coupling between neighboring localized spins oc- 
curs mainly through superexchange processes in the (yz, zx) 
band. The highly anisotropic hopping matrix elements be- 
tween these orbitals, however, leads to an essential depen- 
dence of the superexchange interaction on the orbital con- 
figuration of the minority spin electron (or single hole) of 
each Ru 4+ -ion in the degenerate (xz, yz)-ba,nds. 

In order to gain more insight into the possible form 
of the orbital ordering we performed LDA+U calculations 
for the critical concentration x c = 0.5. Cai.sSro.sRuC^ has 
the space group IA\/acd pq |, i.e. the RuC>6-octahedra are 
only rotated around the c-axis with no tilting (Fig. 0b). 
The RuOg-octahedra remain elongated to the same degree 
as for pure Sr2Ru04 such that the rry-orbital is still higher 
in energy than {xz, yz)-orbitals (Fig. g, right panel). Our 
calculation suggests that, under these conditions, correla- 
tion yields the orbital degeneracy in a (3, 1) state. LDA+U 
calculations generally overestimate the stability of the Mott 
insulating state. In the present case we find for U = 1.5 
eV an insulating groundstate with a charge gap also in the 
(sy)-subband. 

Furthermore, the result of the the LDA+U calculation 
shows that an orbital ordering of the AFO-type is favored 
yielding a ferromagnetic spin exchange. The minority-spin 
electrons occupy alternating xz- and yz-orbitals with a 
slight tilting of the orbital planes away from the c-axis. 
This tilting indicates the partial admixture of d xy -orbitals 
within the LDA+U approach. 



4.2 Region II 

Since magnetic properties are intimately connected to the 
orbital order in {xz,yz}-subband, it is necessary to first 
examine the effect of the crystalline distortion on the or- 
bital order. To this end we performed LDA+U calculation 
for this tilted structure and obtained a rather complicated 
orbital order. The groundstate is an AF insulator. The 
minority-spin electrons (1 per Ru-atom) occupy the or- 
bitals whose planes are in average directed along the a- 
axis (in tetragonal notation (110) direction). However on 
every one of the 4 Ru-atoms in the unit cell those planes 
are rotated from the a-axis by +20° and +15° on one layer 



Fig. 7. Orbital ordering in Cai.sSro^RuCU. At the figure are 
presented minority spin 4-d orbitals of Ru-atoms. z-axis directs 
upwards; x, y-axes is parallel to a bases of oxygen octahedra. 



and by -20° and -15° on the next layer. On the average 
the orbital orientation corresponds to the state d yz + d zx 
Also on one of the 2 Ru-atoms in every layer there is an 
additional tilting of the orbital plane from the long c-axis 
on 34° (see Fig. ??) which corresponds to the admixture 
of the ny-orbital component. The calculation of the easy 
axis using the second order perturbation theory for spin- 
orbit coupling gave the direction of the magnetic moment 
as along the a-axis (tetragonal [110] direction) with a 28° 
tilt from the layer plane (see Appendix). Measurements 
of the uniform magnetic susceptibility show a peak in the 
temperature dependence which is most pronounced for the 
[110] direction of the magnetic field, in agreement with our 
■ LDA+U results M. 



A3 Region I 

The orbital configuration (2,2) is non-degenerate as Hund's 
Rule determines a fully spin polarized S = 1 ionic ground- 
state with a single minority-spin electron in the sy-orbitals. 
The LDA+U method works well here and gives an AF 
Mott insulating groundstate. 

Throughout Region I, there is a first order insulator 
to metal transition as the temperature is raised accom- 
panied by a switch from compressed to elongated RuOg- 
octahedra. It is tempting to interpret this as a first or- 
der transition between the orbital occupancies from (2,2) 
to (3,1). However a detailed examination of this proposal 
has not been made. In this context it would be also in- 
teresting to study the effect of pressure. Since the low- 
temperature phase of Ca2Ru04 has larger volume than 
the high-temperature phase (see Tab. I), pressure tends to 
stabilize the latter. Assuming that indeed the (3,1) elec- 
tronic distribution is realized in this structural state, then 
a orbitally ordered phase accompanied by magnetic order 
based on the localized spin- 1/2 degree of freedom could 
be realized at low enough temperature. Since the crystal 
structure is close to the one at x = 0.2 (see Tab. I) one 
might conclude that AF order would prevail. However, as 
we will see in the next section, the final form of the ex- 
change interactions between the spins is a subtle issue, in 
particular, if a staggered orbital component are included. 



5 Effective Model for the (3,1) electron 
distribution 

5.1 Derivation of spin-isospin Hamiltonian 

In this Section we derive the effective Hamiltonian for the 
localized spin and orbital degrees of freedom in the elec- 
tron configuration (n/ a g\,riy) = (3,1) which covers the 
boundary of region II and III and all of region II. The 
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Table 2. Regions of Ca2- ;E Sr a; Ru04 with orbital occupancy, localized orbital and spin degrees of freedom and order or (dominant 
correlations) (AF = spin antiferromagnetic, FM = spin ferromagnetic, AFO = antiferro-orbital, FO = ferro-orbital). 



Region 


(n y z t zx,n X y) orbital spin order (correlation) 


I (0 < x < 0.2) 

II (0.2 < x < 0.5) 

III (x -> 0.5) 
III (x = 2) 


(2,2) - 5=1 AF 
(3,1) (yz,zx) 5=± FO, (AF) 
(3,1) (yz,zx) S-+1 (AFO) , (FM) 

(1-1) - s=o 



yields antiferromagnetic superexchange through the ex- 
change path, 



nn\-t 

T I 



tin 
o n 



nn 

I T 



This leads to an effective Hamiltonian 



electrons in the d yz - and d za; -orbital are localized, while 
the 7-band remains metallic, although it is half-filled. We 
ignore here this metallic band and concentrate on the lo- 
calized degrees of freedom. The single hole occupying the 
two localized orbitals has a spin-1/2 and orbital index. 
The latter can be described as isospin-1/2 using the no- 
tation |+) for the d yz - and |— ) for the <i 2ar orbital. The 
isospin operator J acts on these states by I z \±) = ±^|±) 
and is a generator of an SU(2) transformation in orbital 
space. 

We will concentrate on the strongest interaction be- 
tween the localized degrees of freedom. The nearest-neighbor 
hopping due to 7r-hybridization between the Ru-d and O- 
p-orbitals lead to the formation of two independent quasi- 
one-dimensional bands, with dispersion in y (x)-direction 
for the d yz (d za; )-orbital. Including the onsite interactions (3) Finally the configuration \+,s) <8> |— ,s') gives an fer- 
the Hamiltonian for these two orbitals has the form romagnetic spin coupling by the following type of path, 



H = J (Si • Si 



-;> 



with an exchange coupling constant 

At 2 
J V 



(5) 



(G) 



(7) 



" ~ * 2-^^ Ci + a x >— ,s Ci >- > s ~^~ C i+a y , + ,x, + ,s C 'h + ,s "•" h"C) 
i,a.s 

+U ^ ni,n,1 n i,n,l + U'J~]ni,+, a ni,-,8' (2) 

i,fi—-\-,— i 

-2J g y^(Sj,+ • Si,- + -n it+ ni-) 
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TTI 
TTI 
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(8) 



Here the intermediate state has lower energy, if the two 
spins are in a triplet configuration due to the Hund's rule 
coupling. The interaction energies for the intermediate 
states corresponding to the spin singlet and triplet state 
are 



where t denotes the hopping matrix element, U and U' arc 
the onsite intra- and interorbital Coulomb repulsion and 
Jh is the Hund's rule coupling (the vector a = (1, 0) and 
(0, 1) connectes the nearest neighbor sites). We can reduce 
the number of parameters by the relation U = U' + 2 Jjj. 
We now consider the case of a nearest-neighbor bond 
along the x-direction and calculation the effective spin- 
spin interaction for different orbital configurations. For 
illustration of the processes involved we introduce the fol- 
lowing self-explanatory notation shown in the two exam- 
ples, 



', T)i ® |+, l)i+a x 
-.T)i® \-,1)i+a !C 



T I 

nn 

t n 
n t 



(3) 

(4) 



The upper (lower) row in the isospinor notation corre- 
spond to the d yz (d Z£C )-orbital. Now we discuss the virtual 
exchange processes leading to the spin interactions for the 
a;-bond. 

(1) The configuration \+,s) <8> |+,s') docs not lead 
to any interaction. (2) The configuration \ — ,s) |— ,s') 



^-Esinglot — U + Jh 
^.Etriplet = U — Jh 



(9) 
(10) 



so that the effective Hamiltonian can be cast into the form 



H" = J'[Si-S 



J i+a x 



1 + a 

2(1 -a) 



with 



./' 



2t 2 J H 



-J 



(3a- l)(a + l) 



(11) 



(12) 



where a = U'/U < 1. Note that a should be sufficiently 
larger than 1/3 in order that the second-order perturba- 
tion approach is valid and the electron distribution (3,1) 
is unique. The orbital configurations in the basis of |+) 
and |— ) are not affected in this second order perturba- 
tion scheme. Therefore the effective isospin Hamiltonian 
is Ising-likc. 

Turning to the bond along the y-axis we need only 
to exchange the role of the two orbitals (d yz and d zx ) 
and obtain the same result. Therefore we can write the 
complete effective Hamiltonian for the localized degrees 
of freedom as, 



V.I. Anisimov et al.: Mott-insulator transition in Ca2- a; Sr :E Ru04: A study of the electronic and magnetic properties 



H = J 



[{A(/? +a + »,.)(/; 



'/a 



Bj S^ +a • Sj 



where the coefficients are given by 
3a 2 + 1 



.4 



B = 



C 



D 



(3a- l)(a + l) 

-(1-a) 2 
(3a 2 + l)(3a-l)(a + l) 
5 — 3a 



4(3a - 1) 



(5-3a)(3a-l) 



»?a = 



?7a = 



(3a-l)(a + l) 2 



2(3a 2 + 1) 



(«x - ««) 



3a 



3a 



/ 2 2 \ 



(13) 

(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



Note that the coefficients rj a and rj' a have opposite sign for 
the x- and y-axis bonds. 

Lattice distortions breaking the tetragonal symmetry 
yield a bias for the local orbital configuration. For the two 
basic orthorhombic distortions, described by the lattice 
strain combinations t\ — e xx —e yy and e 2 = txy (tetragonal 
notation), we can add the following coupling terms to the 
Hamiltonian, 



H dist = ^[ifiei/f + K 2 e 2 If ] 



(20) 




Fig. 8. Mean field phase diagram temperature versus a 
U'/U for x = 0.5. 



effective Hamiltonian for the spin and isospin with mean 
values, 



mi = {If) and s = (S*) 



(21) 



where mi is staggered with opposite sign on the two sub- 
lattices and s is uniform. The coupled self-consistent equa- 
tions are obtained readily, 



mi 



Er=±i re 0Crmi cosh ((3 (A {rg - 



2 , 



B)s) 



2E„ 



: ±1' 



ofiCrmj 



cosh((3(A(7 1 l^^)+B)s 



-22) 



where K\ and K 2 are coupling constants. Note that the 
first term corresponds to a uniform field parallel to the 
z-axis of the isospin, while the second is a transverse field. 
Both drive a ferro-orbital correlation. 



5.2 Mean field discussion 

We analyze now the properties of our effective model by 
means of mean field theory ignoring fluctuations and the 
influence of the itinerant xy-ha,nd. This discussion shows 
that the main features of the different phases where the 
electron distribution (3,1) is applicable are indeed described 
well within our reduced effective model. A discussion be- 
yond the mean field level will be given elsewhere. 



5.2.1 Tetragonal system 

The analysis of the parameters of the effective Hamilto- 
nian shows that a higher energy scale is associated with 
the isospins than with the spins. Their basic interaction is 
of Ising-AFO-type, because C > for all a in the proper 
range. As a consequence the spins correlate ferromagneti- 
cally at a lower energy scale. Hence, we will decouple the 



£r=±i z 0Crmi sinh (g (A (vl - r -^)+B) s ] 
2 Er=±i e &Grmi cosh {p (A (t? 2 



rrni 
2 , 



B)s 



-23) 



where (3 — J/ksT is the inverse temperature in units of 
the exchange energy J. The isospin order has the higher 
transition temperature T determined by the equation 



k B T 



CJ 

4 



(24) 



The onset of FM order occurs at much lower temperature 
so that it is justified to consider mi as already saturated, 
mi — 1/2. Then Curie-temperature for the spin is given 

by 



k B To = —r 



a Tf, 



B 



(25) 



In Fig. @ we show the phase diagram of temperature versus 
the parameter a = U'/U which determines the coupling 
constants. The tendency towards FM order is in good 
agreement with the LDA+U result for x — 0.5, apart 
from the fact that the contributions of a;y-orbital have 
been ignored. This analysis suggests that the orbital or- 
der is governed by the exchange processes included in our 
Hamiltonian. 
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5.2.2 Crystal distortion 

Let us now consider the influence of a spontaneous crystal 
deformation which characterizes region II (0.2 < x < 0.5) 
We will not discuss the origin of the distortion, but take 
it as given with a strength which depends on x. The FO 
bias introduced by a uniform lattice distortion is in com- 
petition with the AFO exchange interaction. In particu- 
lar, a finite strain e2 would introduce a transverse field 
in the effective Ising Hamiltonian for orbital interactions. 
An Ising-model in a transverse field has a quantum critical 
point beyond which the AFO order would be suppressed 
for any temperature. 

We consider here the general case where both ei and 
€2 are finite, as realized in Region II. We distinguish now 
between the two sublattices A and B on the square lattice 
and assign different mean fields to them. Analogous to the 
case above we may separate the spin and the isospin prob- 
lem. First we consider the isospin problem in the absence 
of spin correlation. The corresponding self-consistent equa- 
tions have the form, 



mi a = /(/3,to/ B ) and m IB = f{(i,m IA ) (26) 



with 



/(/3,m) = - 



2Cm + K 1 e 1 



2 v / (2Cm + K^) 2 + K'iej 
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It is easy to verify from this set of equations that for fi- 
nite crystal distortion there is a uniform component m = 
{mi a + tuib)/2 at all temperatures and that the stag- 
gered moment To/ = (mi a ~ "t-/b)/2 occurs only at low 
temperature or not at all, depending on which side of the 
quantum critical point the system lies. Hence the stag- 
gered moment may play a minor role on the background 
of the uniform FO orbital arrangement. 

The numerical solution of Eq. ( J2q , p7| ) is shown in Fig. 9. 
assuming that both e\ and 62 are finite with a ratio K\t\jKi 
3/2. Turning £12 gradually on we observe that the stag- 
gered AFO component is diminished in favor of the FO 
arrangement (Fig. 9a) with a critical point where the stag- 
gered component completely vanishes. Fig. 9b shows the 
orbital orientation on the two sublattices \d) a,b = cos 4>a,b 
sm<f>A,B\—)- Obviously, very close to the quantum crit- 
ical point we find a dominant FO arrangement with a 
weak staggered component. This describes qualitatively 
the favored ordering pattern in the LDA+U calculation 
for x = 0.2, where in each layer the orbital plane has an 
alternating orientation of 20° and 15° with respect to the 
[110] (FO) direction. The comparison with our mean field 
data suggests that the system is then rather close to the 
quantum critical point. 

Assuming again that at the energy scale where spin 
correlation sets in the orbital pattern is basically estab- 
lished, we draw the diagram of transition temperature 



Fig. 9. Effect of lattice deformation K\e\ = 1.5K262 = Ke 
on the orbital order and the effective spin exchange, a) Stag- 
gered orbital moment mi and uniform orbital moment m; b) 
Orientation of the orbital plane cos <f>\+) +sin 4>\~) f° r the sub- 
lattices A and B; c) effective spin exchange ferromagnetic and 
antiferromagnetic 



versus the strain in Fig. 9c. Then we analyze the spin cor- 
relation fixing the orbital order by the T = 0-solution of 
Eq.(E6 27). This is again justified by the observation that 
the isospin mean field has saturated at the temperature 
when spin order sets in. We find weak ferromagnetic cou- 
pling for small distortion and stronger antiferromagnetic 
coupling for larger distortion. The crossover in the spin 
interaction is still in the region of rather strong staggered 
orbital correlations. Note that the crossover point between 
FM and AF spin exchange depends on a and move towards 
larger Ke for smaller a. 

Guided by this phase diagram we may interpret the 
region II in terms of increasing strain as we progress from 
x=0.5 to 0.2. At x=0.2 the quantum critical point has been 
ejipproached but not passed. This interpretation is entirely 
based on the localized orbitals of the a-/3-bands. From 
the LDA+U calculation we conclude, however, that the 7 
-band is also involved in the spin correlations, as is seen in 
the tilting of the orbital plane away from the z-axis for the 
^iifnority spins. As mentioned earlier the 7-band should 
introduce an RKKY-interaction via the Hund's rule cou- 
pling, which is most likely antiferromagnetic, thus further 
diminishing the FM correlation. 

An interesting feature occurs in our effective model 
in connection with the sign change of the spin exchange 
interaction as mentioned above. Ca2_ a: Sr a ;Ru04 as a ran- 
dom alloy does not provide a uniform crystal deformation, 
but the strain depends on the local Sr and Ca distribu- 
tion around each Ru-ion. Thus, in the regime of switching 
between FM and AF coupling, this can lead to frustra- 
tion effects at low temperature. Recent experimental data 
indicate indeed a glass-like behavior [BCJ|. 
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6 Discussion 



A Effect of spin-orbit coupling 



We have presented a consistent picture for the unusual 
phases of the isoelectronic alloy series Ca2^xSr 2 Ru04 based 
on the full multi-band electronic structure (see Table ||). 
Starting from the good metal Sr2Ru04, we find the ef- 
fect of Ca-substitution is to transfer electrons from the 
wider xy-band to the narrower [xz,yz] bands until at a 
critical value of x c — 0.5 there is integer occupancy of 
both subbands. The progressive rotation of the RuC>6 oc- 
tahedra in this region leads to Mott localization of the 3 
electrons in the narrower (xz,yz) bands while the wider 
xy-band which is now half-filled, remains metallic. This 
partial localization of the 4d electrons can explain the puz- 
zling observation of the coexistence of free S = 1/2 local 
moments and metallic behavior in Cai.sSro.sRuO^ The 
actual metal to insulator transition occurs between the 
Ca-rich Regions I (0 < x < 0.2) and II (0.2 < x < 0.5). 
We interpret this transition as driven by a change in the 
orbital occupancy to a completely filled xy-band leaving 
2 electrons in the localized (xz,yz) bands. In this case 
the ordered local spin has the conventional value, S = 1, 
whereas in the metallic Region II short-range correlations 
of a local 5 = 1/2 spin combined with orbital order are 
realized, generating a pronounced anisotropy in the mag- 
netic response. 

One of the curious aspects of Region II is the fact that 
the half- filled xy-ovhital remains itinerant. It is not easy 
to test experimentally which among the bands is really 
responsible for the metallic behavior. The basic feature 
of the state we propose is the presence of a single Fermi 
surface of electron-like character in the xy- or 7-sheet. 
One method which might be able to probe this proposal 
is angle-resolved photo emission spectroscopy (ARPES). 

Recently the electronic structure and magnetic prop- 
erties of Ca2-a;Sr2Ru04 have also been investigated by 
LSDA calculations [£9| . The authors have concluded that 
the crystal structure distortions (rotation, tilt and com- 
pression of the O-octahedra) when Ca is substituted for Sr 
lead to a narrowing and shifting of the bands and hence to 
the enhancement of ferro- and antiferromagnetic instabil- 
ities. However these calculations do not include the effects 
of onsite Coulomb interactions between d-electrons which 
are responsible for the suppression of charge fluctuations 
and for the Mott metal-insulator transition as Ca2Ru04 
is approached. 
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For the Ru 4d-i2g-orbitals spin-orbit coupling plays an im- 
portant role. If we restrict to the three t2g-orbitals, the 
onsite spin-orbit coupling Hamiltonian can is given by 



i a,b,c s,s' 



j 



iCib,, 



(28) 



where A is the coupling constant and e a b c the completely 
antisymmetric tensor (a, b, c — {x,y, z}). We identify the 
indices x, y and z with the orbitals d yz , d zx and d xy , 
respectively. 

In the framework of LDA+U calculation the easy axis 
direction was calculated via minimizing the energy of spin- 
orbit coupling in the second order perturbation theory as 
a function of Euler angles (a, j3, 7) defining the direction 
of magnetization: 
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/ncrk Jn'a'Vi 



TKJ.n (7 
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U^(a,(3,-f)i 



(30) 



where E na (k) is the energy of n-th band for spin pro- 
jection a at point k in the Brillouin zone, / no -k is equal 1 
if the corresponding band is occupied (its energy is below 
the Fermi level) and is equal in the opposite case. &%'<?' 
is the Bloch wave function, l^ and s v are components of 
vector operators of orbital and spin moments, U^ u (a, /3, 7) 
is transformation matrix for rotation of a vector for Euler 
angles (a, £,7). 
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